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^ ■ Abstract 

Location-dependence of physical parameters such as the electromagnetic fine structure constant 
and Newton's G produce body accelerations which violate universality of free fall rates testable 
^vq ' with laboratory and space experiments. Theoretically related cosmological time variation of these 

same parameters are also constrained by experiments such as lunar laser ranging, and these time 
' variations produce accelerations of bodies relative to a preferred cosmological inertial frame. 

cr 

5JT), Experiments which confirm universality of gravitational free fall rates (UFF) also make impressive state- 

ments about the constancy of physical law's key parameters, such as the electromagnetic fine structure 
'constant' and Newton's G. If the mass-energy of a body depends on some parameter a which varies in 
■ space, for example, there is a force on that body 

s dM(a) 2 -. 

b = — c Vina (1) 

oln a 

and a corresponding body-dependent acceleration 63 = F/M which will produce violations of UFF. Objects 
i and j composed of different materials will fall at the differing rates 



dln(Mi/Mj) 2 \Vlna 



dlr 



c2 (2) 



Based on UFF, Einstein made his generalizing hypothesis — the Equivalence Principle — that all physical 
phenomena in local gravity should be identical as that in an accelerating laboratory. His prediction followed 
that two identically constructed clocks located at different altitudes in local gravity must tick at rates dr 
which differ relative to each other as 

d_n = 1 _ a.(n-r 2 ) 

dT2 c z 
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which relationship integrates in a more global Newtonian gravitational potential to be 



dn = 1 U(n)-U(r 2 ) 

dr 2 c 2 

For oscillator atomic clocks whose frequencies are determined by the difference of two quantum energy levels 
of a system, lu^i — (e^ — &i)/h, a location-dependence of a physical parameter a which participates in the 
determination of the energy levels will then affect the Equivalence Principle's prediction of clock rate changes. 
The clock rate relationship given in Equation (4) is altered to now include an additional clock-dependent 
factor 

( dr{\ ( dn\\ ( dln(e kl ) din a \ 

\dr 2 ) kl ~ \dr 2 ) Ep \ dlna dU ) U 

The factor dlna/dU in Equation (5) and the factor \Vlna\/g in Equation (2) are essentially identical, so 
these two equations can be used to compare the ability of clock experiments and UFF experiments to test 
for variation of the parameter a. UFF experiments have been superior performers. Atomic clocks based on 
different hyperfine transitions will typically have factors dln(eki) / din a differing of order one, while different 
elements used in UFF experiments will have nuclear electrostatic energy contributions which render the 
factor din (Mi/Mj)/dln a of order 10 -3 . Clock experiments must therefore be able to measure fractional 
frequency shifts between clocks to precision of a part in 10 10 of the actual gravitational shift in order to 
compete with UFF experiments now performed to the part in 10 13 precision. But the best clock experiment 
has reached only the part in 10 4 precision in testing the EP prediction jlj. Sending atomic clocks on a 
spacecraft which travels to within a few solar radii of the Sun where the gravitational potential grows to 
f CP 6 c 2 could be a competitive experiment if the relative frequencies of different on-board clocks could be 
measured to a part in 10 16 precision. There is no present mission for such an experiment; but some are 
considering it for the future. 



Laboratory experiments which compare the gravitational and inertial forces on differently composed 
material samples, and the ongoing 30 years of lunar laser ranging which compares the free fall rate of Earth 
(f /3 iron-rich core and 2/3 silicate mantle) and differently composed Moon toward the Sun Q are today both 
confirming that gravitational free fall rates are universal to one or two parts in 10 13 precision. A satellite test 
of UFF (the STEP experiment) is under design and development in the United States and Europe to reach a 
part in 10 18 precision, while a quicker, but lower precision, version of such a space mission is being pursued by 
France (MICROSCOPE). UFF experiments are being pushed to these incredible levels of precision because 
no other types of experiments presently compete in their ability to detect location-dependence of fundamental 
physical parameters, or more generally, the presence of any new long range and extremely weak force fields 
in physical law which by coupling to novel attributes of matter produces violations of UFF. 



Lunar laser ranging also places the best present limits on space-time variations of Newton's gravitational 
parameter. If G varies cosmologically in time, the Moon's orbit will show a change in its radius which grows 
linear in time, and phase changes in its various oscillatory motions which grow quadratic in time. Searching 
for the former effect best constrained G in the earlier years of LLR, while now the latter effects dominate. 
G/G < 1CU 12 y^ 1 is LLR's present constraint. 

Theories which predict a non-zero G also generally show that G will be (spatially) dependent on proximity 
to local matter (which produces Newtonian potential U(r)) 

G(r) = Goo fl + V ^r) (6) 
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The coefficient rj vanishes in general relativity but generally is non-zero in alternative metric theories of 
gravity. A body such as Earth with significant internal gravitational binding energy U se i / then experiences 
an anomalous acceleration || || 



dlnM(G) a - / U self 

5a = -^^c VlnG * [l + V^ 1 



(7) 



LLR constrains 77 to be less than 3 10 4 which reflects the fractional precision with which general relativity's 
post-Newtonian structure is confirmed. 

In order to relate location-dependence and cosmic time dependence of a physical parameter, it is consid- 
ered as a function of an underlying dynamical scalar field, a — ► a((f>), with this field having a lagrangian- 
based field equation 

BT 

g^ ^ = 4ttG — = AnG a(4>, r, t) (8) 

For application to experiments looking for violation of UFF on one hand, and to effects related to a cos- 
mological time variation of the physical parameter on the other hand, the scalar field equation takes the 
respective forms 

V 2 (j) = -4irGa{(j), r)/c 2 (9) 
4> +3H<j) = AkG a((f>, t) (10) 

with H = R/ R being the Hubble expansion rate of the universe, R(t) being the size parameter of the 
universe. A quantitative relationship can then be sought between violations of UFF and time evolution 
of the physical parameter. For two materials i and j the anomalous forces given by Equation (2) lead to 
acceleration differences 

\di — 3j\ dlna dln(Mi/Mj) 



\g s \ d<j) dh 



a s dV j \ p s dV (11) 



with p s and a s being the mass density and scalar density, respectively, of the source(s) which accelerates 
bodies i and j of respective mass-enrgies Mj and Mj. On the other hand, a first integration of Equation 
(17) yields an expression for the cosmological time evolution of a 

la f 1 a\ ( R 3 H \ dlna/d<f> 



Ha \H aj po \dlna/d(/)J pg R 3 H 

+ L5 ^y*vy ^% dpl (12) 

d<t> J Po R 3 H p(p>) c + (A + cf>2)/3irG 

with subscript c indicating averaged cosmological densities of the scalar and mass-energy sources, and A 
being the cosmological constant. A convenient variable dp = dRj R = H dt has been used to parameterize 
cosmological time evolution, p Q is a useful historic epoch and p is the present (unlabeled quantities are their 
present values). 

As example, consider the fine structure parameter a = e 2 /hc. As previously discussed, UFF has been 
confirmed to about a part in 10 13 among elements whose nuclear electrostatic energy contributions to mass 
energies typically vary by parts in 10 3 of the whole 

0'"(AW) „ 10 _ 3 (13) 
ulna 
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Equation (11) can then be used to eliminate the numerically unknown factor dlna/d<j) from the cosmological 
time evolution Equation (12) to produce the constraint 

I* < 10-™ f p ( R 3 (p')h(p') \ ( p(p) s v(p')c \ dp' 

Ha ~ J Po \ R 3 (p)H(p) J \a(p) s p(p')J \ + {k + &)/^Gp(p') 

in which p a is here chosen to be a past epoch when a was zero. No other direct observations have constrained 
a/a to near the level of 10 -10 H. This raises the question of whether circumstances exist for which the 
cosmological time evolution observations produce interesting constraints on physical theory? The first factor 
in the integrand is less than one in our expanding universe. If the scalar density is an attribute of ordinary 
matter and concentrates in bodies along with ordinary matter, the second factor should also be of order 
one unless there is something very unusual about the present epoch. In the "attractor" scenario for scalar 
theories || 0], it is found that the dynamical, cosmological background scalar field will evolve toward a 
value where it's coupling function to its source, /(</>), has a minimum (and hence zero derivative) — if such 
locations exist in /(</>). If at the present epoch this field is very close to this site, then a(p ')/a(p) could be 
a huge factor which compensates for some or much of the 10 -10 factor in Equation (14). Another way in 
which the integral in this equation could have a very large value is if the cosmological scalar source has little 
or nothing to do with ordinary matter and does not concentrate in ordinary matter sources such as stars or 
planets to the degree that the ordinary matter does . 



Cosmological time variation of a physical parameter also establishes a cosmic preferred frame and body- 
dependent Aristotelian-like accelerations ||. If one performs a Lorentz transformation to the rest frame of 
a body which moves at velocity v relative to the locality's special frame in which a parameter has pure time 
dependence, a spatial gradient is acquired in the new frame, a(t) — > a(t' + v ■ r'/c 2 ). From equation (2), a 
force then accelerates the body in proportion to its velocity in the cosmos 

dv a dlnM _, 

7i7 ~ai v ( 15 -> 

at a aina 

Another way to understand the origin of this acceleration rests on the time variation of the body's mass and 
conservation of momentum 

= !L(M(t)v) = M'v + Mv = Mtj + --^-v (16) 
dt a ulna 

which reproduces Equation (15). 



This bi-product of a cosmically time- varying physical parameter alters the response of a binary pulsar 
system's orbit to a changing gravitational coupling strength G(t). A near-circular binary system has period 
P = 2ttL 3 /(GM) 2 plus relativistic corrections, L is the angular motion quantity \ f x v\, and M is the mass 
of the system. Under changing G both L and M change as well, so 

P L G M 

P^ 3 L- 2 G- 2 M < 17 > 

For gravitationally compact bodies (those having large dlnM/dlnG) such as neutron stars the fractional 
time rate of change of both I and m can be a few tenths of that for G and have signs working in opposition 
to G. The P/G ratio is therefore generally significantly reduced for such systems. 



If the two bodies in the binary system have different gravitational compactnesses, then the system will 
also experience a polarizing force in the direction of the system's velocity relative to the cosmic preferred 
frame, and the resulting orbital perturbation can also be sought as a measure of time variation of G. 
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